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We show that the right-handed (RH) sneutrino in the NMSSM can account for the observed excess 
in the Fermi-LAT spectrum of gamma-rays from the Galactic Centre, while fulfilling all the current 
experimental constraints from the LHC as well as from direct and indirect dark matter searches. 
We have explored the parameter space of this scenario, computed the gamma-ray spectrum for 
each phenomenologically viable solution and then performed a y 1 2 fit to the excess. Unlike previous 
studies based on model-independent interpretations, we have taken into account the full annihilation 
spectrum, without assuming pure annihilation channels. Furthermore, we have incorporated limits 
from direct detection experiments, LHC bounds and also the constraints from Fermi-LAT on dwarf 
spheroidal galaxies (dSphs) and gamma-ray spectral lines. In addition, we have estimated the effect 
of the most recent Fermi-LAT reprocessed data (Pass 8). In general, we obtain good fits to the 
GCE when the RH sneutrino annihilates mainly into pairs of light singlet-like scalar or pseudoscalar 
Higgs bosons that subsequently decay in flight, producing four-body final states and spectral features 
that improve the goodness of the fit at large energies. The best fit (y 2 = 20.8) corresponds to a 
RH sneutrino with a mass of 64 GeV which annihilates preferentially into a pair of light singlet¬ 
like pseudoscalar Higgs bosons (with masses of order 60 GeV). Besides, we have analysed other 
channels that also provide good fits to the excess. Finally, we discuss the implications for direct 
and indirect detection searches paying special attention to the possible appearance of gamma-ray 
spectral features in near future Fermi-LAT analyses, as well as deviations from the SM-like Higgs 
properties at the LHC. Remarkably, many of the scenarios that fit the GCE can also be probed by 
these other complementary techniques. 


INTRODUCTION 


A large number of cosmological and astrophysical ob¬ 
servations have evidenced that 85% of the matter con¬ 
tent of the Universe is in the form of dark matter (DM). 
A generic weakly interacting massive particle (WIMP) 
is a well-motivated candidate for this new kind of mat¬ 
ter, since its thermal production in the Early Universe 
would match the observed DM abundance. In addition, 
WIMPs can be easily accommodated in theories beyond 
the Standard Model, such as Supersymmetry (SUSY). 

WIMPs can be searched for indirectly, through the 
particles produced when they annihilate in the DM halo 
(photons, neutrinos and antiparticles). Among the dif¬ 
ferent annihilation products, gamma-rays provide an ap¬ 
pealing detection possibility because the signal can be 
traced back to the source. The Large Area Telescope 
(LAT) aboard the Fermi Gamma-ray Space Telescope has 
produced the most detailed maps of the gamma-ray sky 
for a wide range of energies, with unprecedented angular 
and energy resolutions. Using data from the Fermi-LAT, 
various studies have revealed the presence of an excess 
from an extendedgamma-ray source in the inner region 
of the Galaxy [lHfJ, a signal that is robust when known 


uncertainties are taken into account EM- Although 
the explanation of this Galactic Centre excess (GCE) is 
still under debate], if it were interpreted in terms of DM 
annihilations dl 1T3| it would correspond to a particle 
with a mass in the range 30 — 50 GeV for a bb final state 
(7—10 GeV for a rf final state) and with an annihilation 
cross section in the DM halo, {<jv)o ~ 1 — 2 x 10 -26 cm 3 /s, 
remarkably close to that expected from a thermal relic. 

Several attempts have been made to explain the GCE 
in terms of simplified models for DM 12|, ll4H5l| . con¬ 
sidering DM annihilation into pure channels. However, 
as pointed out in Ref. |52j], it is crucial to investigate if 
this excess can be obtained within a complete theoretical 
framework. For example, it has been recently shown that 
the neutralino could reproduce the GCE for DM masses 
up to hundreds of GeV depending on the primary annihi¬ 
lation channel within the context of the MSSM 11] and 
the Next-to-MSSM (NMSSM) [5^®. 

In this article, we carry out a complete analysis of the 
right-handed (RH) sneutrino in the NMSSM [53, @| and 


1 A recent study shows that the GCE might be produced, within 

current uncertainties, by a population of unresolved millisecond 
pulsars flU . 
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demonstrate that it can successfully account for the GCE 
while fulfilling all the experimental constraints from di¬ 
rect and indirect DM searches as well as collider physics. 
We apply the LUX and SuperCDMS limits on the spin- 
independent elastic scattering cross section of DM off 
protons, which are currently the most stringent bounds 
from direct detection experiments. We also consider the 
latest results from the LHC on the Higgs boson mass and 
couplings to the SM particles, which are known to be spe¬ 
cially constraining for light DM scenarios through the 
upper bound on the invisible and non Standard Model 
Higgs decays. Besides, the latest bounds from the mea¬ 
surement of the rare decays B$ —*■ y + , b —> s'y and 
B+ T+ v T are also applied. Finally, we incorporate 
the Fermi-LAT constraints on dwarf spheroidal galaxies 
(dSphs) and spectral feature searches in the gamma-ray 
spectrum, including an estimation of the effect that the 
most recent results derived from the Pass 8 data impose 
on our results. 


FITTING THE GCE WITH RH SNEUTRINOS 


This model has been extensively described in Refs. [57, 
[Hit ]. It is an extended version of the NMSSM, in which 
a new gauge singlet superfield, N, is introduced in order 
to account for RH neutrino and sneutrino states as in 
mm. The superpotential of this construction is given 

by 


W = U'nmssm + XnSNN + y N L ■ H 2 N , 


( 1 ) 


where flavour indices are omitted and the dot denotes the 
antisymmetric SU (2) l product. Wnmssm is the NMSSM 
superpotential, Xn is a new dimensionless coupling, yw 
is the neutrino Yukawa coupling, and L/ 1,2 are the down 
and up type doublet Higgs components, respectively. As 
in the NMSSM, a global Z 3 symmetry is imposed so that 
there are no supersymmetric mass terms in the superpo¬ 
tential. Since we assume 17-parity conservation in order 
to guarantee the stability of the LSP, the terms NNN 
and SSN are forbidden. Furthermore, we do not con¬ 
sider CP violation in the Higgs sector. 

After radiative electroweak symmetry breaking the 
Higgs fields get non-vanishing vacuum expectation val¬ 
ues (VEVs) and the physical Higgs states correspond to a 
superposition of the Hd, H u and S fields. The RH sneu¬ 
trino interacts with the SM particles through the mix¬ 
ing in the Higgs sector thanks to the coupling XnSNN, 
thereby behaving as a WIMP. 

Interestingly, light RH sneutrinos with masses in the 
range of 10 — 150 GeV are viable DM particles [Hl| and 
constitute ideal candidates to account for the GCE, as 


we already pointed out in Ref. [50]. Their phenomenol¬ 
ogy is very rich, as they can annihilate into a variety 
of final states, some of which include scalar and pseu¬ 
doscalar Higgses. In particular, if > m H o^cq, the 


annihilation final state of RH sneutrinos is dominated 
by a HiHi(AiAi) pair in vast regions of the parameter 
space. The subsequent decay of each scalar and pseu¬ 
doscalar Higgs into pairs of fermions or gauge bosons 
gives rise to non-standard final states, which often dis¬ 
play spectral features coming from the 77 final states. 

Given that the RH sneutrino annihilation contains a 
mixture of final states, often including exotic configura¬ 
tions, the GCE model-independent approach generally 
found in the literature is not applicable To fit the GCE 
we have followed the approach of Ref. [Toj , where the au¬ 
thors take into account theoretical model systematics by 
exploring a large range of Galactic diffuse emission mod- 
el4j. Regarding the DM distribution, we have considered 
a generalised NFW profile with slope 7 = 1.2 in the re¬ 
gion of interest (ROI) 2° < | 6 | < 20° and \l\ < 20°, as in 
Ref. 0. 

To implement the aforementioned analysis, we have 
performed a series of scans over the parameter space of 
the model in the RH sneutrino mass range 1 — 150 GeV, 
computing the gamma-ray spectrum as well as the RH 
sneutrino relic abundance with micrOMEGAs 3.6.9 163 ] . 
The number of free parameters and the correspon ding 
ranges of variation coincide with those used in Ref. [56| 
(only the range in the RH sneutrino mass has been en¬ 
larged to accommodate new solutions). In our anal¬ 
ysis, we have also incorporated the most recent con¬ 
straints from direct detection experiments using the lat¬ 
est results of LUX and SuperCDMS 63 65|. LHC con¬ 
straints on the masses of supersymmetric particles, the 
mass and couplings of the SM Higgs boson, as well as 
bounds from the rare decays B$ — > y + , b —> s'y and 
B + —> r + py, have also been implemented (for more de¬ 
tails, see Ref. [5(1). 

We have set an upper bound on the RH sneutrino 
relic abundance, Djo/i 2 < 0.13, consistent with the lat¬ 
est Planck results [ 66 |. Besides, we have considered the 
possibility that RH sneutrinos only contribute to a frac¬ 
tion of the total relic density and set for concreteness 
a lower bound on the relic abundance, 0.001 < ^JVi h 2 - 
To deal with these cases, the fractional density, £ = 
min[l,D J 7 i /i 2 / 0 . 11 ], has been introduced to account for 
the reduction in the rates for direct and indirect searches 
(assuming that the RH sneutrino is present in the DM 
halo in the same proportion as in the Universe). 

Then, we have convoluted the differential photon spec¬ 
trum of the scan points fulfilling all these constraints with 


2 When these model systematics are included as correlated errors 
in the residual spectrum, the best fit for the DM interpretation 
of the model-independent analysis of Ref. is obtained for a 
bb final state with a mass of 49+54 GeV and a velocity averaged 
annihilation cross section of 1 . 76 ^q 27 x 10 —26 cm 3 /s. Other 
analyses of the GCE employ different assumptions on the Galac¬ 

tic diffuse and point source components, and the reconstructed 
DM mass and annihilation cross section differ slightly. 
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FIG. 1. RH sneutrino annihilation diagrams that produce lines and box-shaped features in the gamma-ray spectrum. 


the energy resolution of the LAT instrument. We have 
used the P7REP-S0URCE-V15 total (front and back) reso¬ 
lution of the reconstructed incoming photon energy as a 
function of the energy for normally incident photons. 

Afterwards, we have calculated the y 2 function as fol¬ 
lows [loj : 



where £y is the covariance matrix containing the statis¬ 
tical errors and the diffuse model and residual system- 
atics 0- dN/dEi ( dN/dEi ) stands for the measured 
(predicted) flux in the ith energy bin. The vector 6 refers 
to all parameters of our model which determine the pre¬ 
dicted flux. 


Constraints from indirect DM searches 

The Fermi-LAT satellite has also provided bounds 
on the DM annihilation cross section in the Galactic 
halo derived from the study of the gamma-ray spectrum 
from dSphs and the search for spectral features in the 
Galactic Centre. These limits play an important role in 
the current analysis. Let us review them in more detail. 


• Dwarf spheroidal galaxies 

The mass of these objects is dominated by DM; hence, 
they constitute ideal targets for indirect searches. The 
Fermi-LAT Collaboration has performed an analysis of 
the gamma-ray emission from 25 dSphs using four years 
of data [67j. The absence of a signal can be interpreted 
as constraints on the annihilation cross section of DM 
particles. It is customary to assume annihilation into 
pure SM channels in the calculation of these bounds. 

The occurrence of non standard annihilation final 
states in our model prevents us from using these results 
directly. Instead, we have extracted independent upper 
bounds on £, 2 (crv}o for each of the six more constraining 
dSphs (Coma Berenices, Draco, Segue I, Ursa Major II, 
Ursa Minor and Willman I), using the DM flux predicted 
by our model and the mean values for the J factors from 


Ref. [67J. Then, we have applied the most restrictive 
of these limits to our data. We have checked that this 
method leads to slightly less stringent bounds than the 
combined limit from the Fermi-LAT Collaboration (by a 
factor smaller than 1.5), when applied to the region of 
DM masses from 10 to 100 GeV with pure annihilation 
channels. 

Lastly, we have also estimated the impact of the 
preliminary results derived from the latest data (Pass 8) 
presented by the Fermi-LAT Collaboration [68[. In 
general for any final state, the limit on £ 2 (ctt;}o improves 
by approximately a factor 4 — 5 for a DM mass in the 
range 10 — 130 GeV. Conservatively, we have used a 
factor 4 to assess the dSph bounds derived from the 
newest Fermi data. 


• Gamma-ray spectral features 

The Fermi-LAT Collaboration has performed a search 
for spectral lines in the energy range 5 — 300 GeV [69| . 
Not having found any globally significant spectral fea¬ 
ture, this analysis has been translated into 95% C.L. up¬ 
per limits on the DM annihilation cross section into a 
pair of photons, (cru)o, 77 . 

The RH sneutrino in the NMSSM can give rise to a 
complex spectrum, displaying lines and box-shaped spec¬ 
tral feature^, which arise from the diagrams shown in 
Figure[T] and involve annihilation into pairs of light scalar 
and pseudoscalar Higgs bosons. The first diagram shows 
the usual contribution to the primary production of a 
pair of photons through a loop of charginos, sfermions, 
top quarks, W ± and charged Higgses, and would pro¬ 
duce a line with E 1 = . The second and third dia¬ 

grams would produce a line with energy U 7 = m H ° fc /2 
or Ery = m A o 2 /2 if any of the Higgs bosons were pro¬ 
duced nearly at rest, when 2 m Ali ss m H o + m H o or 
2 m<j « m ,to + m/io. Otherwise, the decay in flight 
of the boosted Higgs bosons gives rise to box-shaped 
features with a maximum energy E™ ax = m^/2 and 
widths A Ery = 77 — m H o k (m A o 2 ). 


3 For a more detailed discussion on these features see Refs. mm- 
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FIG. 2. Velocity-averaged annihilation cross section of RH sneutrinos in the Galactic halo as a function of the RH sneutrino mass. All the 
points provide a fit to the GCE at 95% C.L. We have also imposed all the experimental bounds, including dSph constraints and Fermi-LAT 
searches for spectral lines. The colour code indicates different dominant final states of the RH sneutrino annihilation. The best fit points 
for each annihilation channel are represented by a star. Points allowed by our estimation of the Pass 8 dSph bounds are encircled in black. 


The published bounds [69] on (av) o )77 do not include 
the specific DM halo used in this paper for the anal¬ 
ysis of the GCE. In order to recalculate this limit, we 
have computed the J factor for our halo in the region of 
interest (ROI) R410, Jgce = 13.0 x 10 22 GeV 2 cm -5 , 
and compared it with the one used in Ref. (69| . which 
yields Jlat = 8.53 x 10 22 GeV 2 cm -5 . The ratio 
r = Jlat/Jgce ~ 1-52 is then applied to the Fermi- 
LAT bounds on (av) o i77 from Ref. |69|. 

Then, we have applied the bounds on the annihilation 
cross section into two photons to the £ 2 (av) o )77 predicted 
by our model for monochromatic gamma-ray lines. Con¬ 
cerning the box-shaped contributions, first we have de¬ 
rived the corresponding limits on the annihilation cross 
section from the Fermi-LAT gamnra-ray line boundfH 
and afterwards we have applied them to our prediction 
weighted by the fractional DM density squared, £ 2 , along 
the box width. Finally, note that stronger bounds could 
be obtained if we used a ROI which is optimised for the 
profile used here, but this is out the scope of this article. 


4 This ROI is defined as a 41° circular region centred on the Galac- 
tic centre with a mask applied to | 6 | < 5° and \l\ > 6 °, and has 
been optimised for a regular (7 = 1 ) NFW profile. 

5 The use of this derived bound for box-shaped features is well- 
motivated since the energy binning of the Fermi-LAT flux is cho¬ 
sen to be of the order of the energy resolution of the instrument. 
Hence, we are allowed to approximate this contribution as a con¬ 
tinuum of lines extending from the minimum to the maximum 
box energy. 


RESULTS 

In Figured we show £ 2 (av)o vs. for the points of 
the parameter space that fit the GCE at 95% C.L. The 
different colours indicate the main annihilation channel 
(remember that the whole annihilation spectrum is con¬ 
sidered when calculating the gamma-ray flux). The stars 
represent the best fit point for each of the dominant an¬ 
nihilation channels and their properties are summarised 
in TableO] where we distinguish “pure final states” (if the 
main annihilation channel contributes to more than 90% 
to £ 2 (<7 u)o) and “mixed final states”. Lastly, black circles 
correspond to the points that would be allowed by our 
estimation of the Pass 8 constraints on dSphs. 

As we can observe, there are solutions that fit 
the GCE for RH sneutrino masses in the range 
rrijq = 15—135 GeV, while fulfilling all other experi¬ 
mental constraints (from direct and indirect dark matter 
searches as well as from the LHC). The best fit points for 
pure annihilation channels are i n g ood agreement with 
model-independent studies Hlf but we have also 
obtained new non-standard annihilation channels (into 
light scalar and pseudoscalar singlet-like Higgs bosons) 
and examples with mixed final states which provide 
a slightly better fit to the GCE. Notice also that the 
points are separated in two regions in the RH sneutrino 
mass. Let us comment in more detail these two regions. 

• mjy ~ 15 — 30 GeV. We have found solutions where 
the RH sneutrino annihilates mainly into a pair of 
very light, singlet-like, CP-odd Higgs bosons (cyan). 
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FIG. 3. Differential gamma-ray spectrum for the points in Table [I] for pure (left panel) and predominant (right panel) annihilation 
channels. The colour convention is as in Fig. [2] The experimental data and errors are extracted from Ref. na , as well as the best fit for a 
pure bb channel, represented by a black dashed line. 


Pure final states 


Final state 

m N 1 (GeV) 

£ 2 (av)o (cm 3 /s) 

^ N ^ 2 

x 2 

H°H° (91.8%) 

119.8 

5.1 x 10~ 26 

0.094 

21.9 

A°A° (90.6%) 

65.0 

2.7 x 10~ 26 

0.109 

22.3 

bb (90.2%) 

46.1 

1.9 x 10 -26 

0.038 

22.6 


Mixed final states 


Final 

state 

m N 1 (GeV) 

t 2 {av) 0 

(cm 3 /s) 

JO 

10 

x 2 

A°Ai 

(44.7%) 

63.8 

2.9 

X 

10“ 

26 

0.061 

20.8 

bb 

(42.1%) 

63.2 

2.9 

X 

10' 

26 

0.042 

21.0 

H°Hi 

(71.4%) 

121.4 

5.4 

X 

10' 

26 

0.075 

21.6 

99 

(38.8%) 

39.6 

1.4 

X 

10' 

26 

0.071 

23.7 

CC 

(33.0%) 

39.0 

1.2 

X 

10' 

26 

0.099 

25.4 

mm 

(44.5%) 

127.4 

4.3 

X 

10' 

26 

0.054 

25.9 

AiA° (4r) 

(67.5%) 

25.5 

1.5 

X 

10' 

26 

0.068 

27.4 

W+W~ 

(28.0%) 

72.4 

2.6 

X 

10' 

26 

0.104 

29.2 


TABLE I. Properties of the points that provide the best fit 
to the GCE for different annihilation final states. We have 
separated the solutions into pure final states (which have an 
annihilation percentage into a given channel bigger than 90%) 
and mixed final states (in which case we show the dominant 
channel with its percentage). 


Since m^o < 2 mb, these pseudoscalars cannot decay 
into a pair of b quarks and instead they do it predomi¬ 
nantly into a pair of r leptons. The resulting process, 


N\Ni —> 2 A® —> 4t, leads to a leptonic final state 
(with best fit around ~ 25 GeV), which differs 

from the usual 2 r final state (whose best fit is around 
10 GeV 0 ). We have also found 2 t final states; 
however, these appear only for < 5 GeV {56] and 
therefore fall out of the 95% C.L. 

• niff ~ 30 — 135 GeV. This region is populated by 
points which present annihilation mainly into bb (grey), 
cc (green), gg (violet), (cyan), (blue) and 

H°H$ (dark blue). 

The best fit for a pure annihilation into a bb pair is 
obtained for m# = 46.1 GeV (see Table [TJ), in good 
agreement with Ref. 0 , but it shifts to larger masses 
m j y i = 63.1 GeV if mixed final states are considered. 
Very few solutions with dominant cc and gg final states 
are found. These channels dominate when the up com¬ 
ponent of the lightest Higgs is larger than the down com¬ 
ponent, which enhances the Higgs coupling to up-type 
fermions and top loop contributions to gg final states. 
However, these loop contributions also enhance the 77 
line production and most of the points are excluded 
for this reason. Besides, these final states are always 
related to the resonant annihilation of RH sneutrinos 
through a light singlet-like H° [56] and typically have 
a smaller relic abundance than the lower bound consid¬ 
ered in this article. This also happens for other channels 
when ~ m H a/2 k, 63 GeV and explains the gap in 
the plot. 

The annihilation into a pair of CP even Higgs bosons 
takes place mostly for >60 GeV. These subse¬ 
quently decay mainly into bb (if the down component 
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of H® is large), giving rise to a four quark final state 
NiNi —> 2 Hi —> 46. It can also decay into gluons and cc 
(if the up component of H° is large). It is worth noting 
that this light H® is singlet-like, typically has a mass of 
30 — 100 GeV, very typical of the NMSSM. For this rea¬ 
son the best fit point (which in our analysis is obtained 
for rn j y i « 121 GeV with y 2 = 21.6) differs from other 
solutions found for the SM Higgs [121 ]. 

Finally, regarding the annihilation into A®A®, the 
pseudoscalar in this region satisfies m^o > 2 mb and thus 
decays predominantly into 66. The annihilation into four 
6 quarks, NiNi —»■ 2 A® —> 46, produces a best fit of the 
GCE for = 65.0 GeV when this channel is almost 
pure. Mixed scenarios are in fact responsible for the best 
fit point found in our analysis with 44% of annihilation 
into A®A®, a RH sneutrino mass m i y. | =63.5 GeV and a 

X 2 = 20.8. 

The spectra for the best fit points in Table |I] are shown 
in Figure [3] As we can see the pseudoscalar and scalar 
channels provide the best fits. This is mostly due to 
the presence of spectral features (mainly lines in the case 
of A® A® and box-shaped emissions for H®H®) that im¬ 
proves the fit at large energies. On top of this, the en¬ 
ergy distribution of the spectra of four-body final states 
(fermions, gauge bosons or combinations between these 
two, as shown in Fig.l) is also known to provide better 
fits to the GCE [37| . Other two-body final states, such 
as 66, cc, gg and W + W ~, also benefit from the pres¬ 
ence of lines in the spectrum which might be probed by 
Fermi-LAT using the latest Pass 8 reprocessed data. 

In Figure EJ we also show all the points that are not 
ruled out by our estimation of the Pass 8 dSph bounds 
and provide a 95% C.L. fit to the GCE. It is remark¬ 
able that all the solutions found that fit the GCE at 68% 
C.L. would be now ruled out, highlighting the tension 
between the Pass 8 data and the GCE. Moreover, the 
points that survive the Pass 8 constraints are mainly re¬ 
lated to dominant annihilation into non standard final 
states, such as A^A® and . This is due to the fact 

that these solutions produce four particles in the final 
state, which makes their spectrum different from that 
of the usual two-body final states. Furthermore, at the 
same time they produce spectral features at high ener¬ 
gies. It is worth noting that the points which are still 
alive and fit the GCE at 95% C.L. might show up in the 
results of the Fermi-LAT Collaboration in the very near 
future. 


Implications for direct DM searches 

The RH sneutrino elastic scattering off quarks is me¬ 
diated by the exchange of a CP-even Higgs boson on 
a t channel, and the dominant contribution arises from 


the lightest of these particles, H°. As it was shown in 
Ref. [61], the annihilation cross section and the spin- 
independent scattering cross section off protons, oA 1 , 
are correlated by crossing symmetry when the main an¬ 
nihilation final state in the Early Universe is any quark 
pair. This generic.ally leads to a large , already ex¬ 
cluded by current direct detection limits. However, this 
correlation can be broken in the presence of resonant an¬ 
nihilations. Because of this as well as the different anni¬ 
hilation final states, the predictions for eA 1 span many 
orders of magnitude [HU. 

As already emphasised in the previous section, there 
is a wide variety of final states that can fit the GCE at 
95% C.L.; therefore the prospects for direct detection can 
be very different. In Figure |4] we show the theoretical 
predictions for as a function of the RH sneutrino 

mass for all the points that fit the GCE at 95% C.L. We 
have also included the most stringent bounds from direct 
detection experiments 63 - 65| and some of the predicted 
sensitivities of next-generation detectors. 

It is noteworthy that many of the points have a rela¬ 
tively large scattering cross section and might be within 
the reach of future experiments, such SuperCDMS and 
LZ. In particular, this is the case of most of the exam¬ 
ples with A^A® final states and >30 GeV that also 
satisfy our estimation of the Pass 8 constraints. On the 
contrary, especially in the points with resonances, the 
predicted £ a ^ ip can be extremely small. The resonance 
with the SM-like Higgs Ul) is narrow and occurs around 
to i y | « m H o/2 ss 62.5 GeV. The resonance with the 
lightest singlet-like Higgs, H°, seems broader due to the 
variation of rn H o throughout the scan and generically oc¬ 
curs for <60 GeV. Fortunately, the presence of such 
light scalar and/or pseudoscalar Higgs bosons in these 
cases might give rise to interesting signals at the LHC. 


Implications for indirect DM searches 


As already pointed out, most of the cases that pro¬ 
vide good fits to the GCE also contain box-shaped fea¬ 
tures and lines in the gamma-ray spectrum. Although 
we have included in our analysis the constraints from 
the Fermi-LAT Collaboration on these contributions, fu¬ 
ture searches for these kind of features might provide 
new insight on the DM properties. In this subsection, we 
will determine which of the points that fit the GCE at 
95% C.L. are also within the reach of future searches for 
gamma-ray lines. 

To this aim, we have defined the following ratio: 


{av)^ T {E 1 )/r 


( 3 ) 


which quantifies how far our predictions are from the 
current experimental limits. Here ^ 2 (av) s f is the contri- 






7 



10 10 2 

rrv (GeV) 


m~ (GeV) 


m a (GeV) 


w w 


FIG. 4. Theoretical predictions for as a function of for points which fit the GCE at 95% C.L. and fulfil all the experimental 

bounds. The colour convention is as in Fig. [2] Solid lines represent the current experimental upper bounds from direct detection experi¬ 
ments, whereas dotted lines are the projected sensitivities of next-generation detectors. The dashed line corresponds to an approximate 
band where neutrino coherent scattering with nuclei will begin to limit the sensitivity of direct detection experiments. Closed contours 
represent the areas compatible with the observed excesses in DAMA/LIBRA (orange), CRESST (red), CDMS II (blue), and CoGeNT 
(green). 
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FIG. 5. R-yy ratio as a function of the photon energy, E-,. for 
which the ratio is maximised. All the points fulfil all the experi¬ 
mental bounds. The colour convention is as in Fig. [2] The dashed 
line denotes our estimation of the improved sensitivity to spectral 
feature searches with Fermi-LAT Pass 8 data. 


bution from our model to spectral features (box or lines) 
on the gamma-ray spectrum and (<t?;} 7 / 4t is the expected 
limit from Ref. [69| for R41. The factor r « 1.52 has been 
introduced in the previous section to convert the bounds 
of Ref. [69} to the DM halo considered in this article. 
i? 77 is evaluated at the energy Ej, for which the ratio is 


maximised. For gamma-ray lines E 1 coincides with the 
energy of the line. On the other hand, for box-shaped 
features E 7 represents the mean value of the Fermi-LAT 
energy bin for which £ 2 (or’) s / is closer to the current 
expected limit. Notice that i? 77 < 1, since the current 
bound on spectral features has already been applied to 
our data. 

In Figure O we represent i? 77 vs. E 1 for all the points 
that fit the GCE at 95% C.L. We also indicate with a 
dashed line the expected improvement on these kind of 
searches with the Pass 8 dated. We can observe that 
many of these scenarios have i? 77 > 0.5, which means 
that an improvement of a factor 2 in the Fermi-LAT 
sensitivity to the search for spectral features would be 
enough to probe these solutions. 

As has already been emphasised, the points in which 
the RH sneutrino annihilates mainly into scalar or pseu¬ 
doscalar Higgs bosons typically present box-shaped fea¬ 
tures and/or lines in their spectrum. The energy E 7 at 
which we evaluate I? 77 (and at which we expect them 
to be detected) is systematically shifted towards low val¬ 
ues, since the sensitivity of the Fermi-LAT instrument 
is betteiQ. In particular, all the points with H®Hi fi¬ 
nal states have E 7 < 60 GeV (in spite of the RH sneu- 


6 This estimation is obtained assuming that the improvement on 
the expected limit with respect to the Pass 7 data for the Einasto 
profile (ROI R16) shown in Ref. [7^1 can also be applied to the 
DM halo considered in this article. 

7 When both lines and box-shaped features are present in the same 



























BR(H°->inv) 


appear as non-standard decay channel^ i7 2 H}H} 
and H} -f A}A}). 

Current LHC data leave considerable room for these 
contributions [e.g., BR(7J 2 —> inv) < 0.27]. Nevertheless, 
the future high luminosity LHC run (3000 fb -1 at y's = 
14 TeV) will significantly narrow down these ranges (e.g., 
the invisible branching ratio in ZH associated production 
will be probed at 95%C.L. down to 13% and the Higgs 
couplings to SM particles will be measured individually 
with 10% accuracy jz3). It is therefore conceivable that 
some of the scenarios considered in this article can be 
probed in this way. 

We have defined the invisible branching ratio of the 
SM Higgs and the branching ratio of the SM Higgs into 
non SM particles as follows, 


FIG. 6. Branching ratio of the SM-Higgs, H}, into non SM par¬ 
ticles as a function of its invisible branching ratio. All the points 
fulfil all the experimental bounds. The colour convention is as in 
Fig-11] The regions above the dashed line and to the right of the 
dot-dashed line will be probed by future searches at the LHC. 


trino mass being ft 60 — 135 GeV). Regarding the 
points with A}A} final states, only a portion of them with 
E-y ft 70 GeV might be observable in the near future. 

Unfortunately, all of the points that satisfy our expec¬ 
tation of the Pass 8 bound on dSphs have very small 
values of i? 77 and therefore would not present any ob¬ 
servable feature in the gamma-ray spectrum. 


Implications for LHC searches 

Finally, in this section we investigate the implications 
of new physics searches at the LHC. In general, the 
NMSSM introduces profound changes in the structure 
of the Higgs sector. Most notably, the addition of new 
states opens the door to the occurrence of light singlet¬ 
like scalar and pseudoscalar Higgs bosons with a mass 
that can be considerably smaller than that of the SM 
Higgs. As we have seen in the previous sections, these 
scenarios are very common when we try to fit the GCE 
with RH sneutrino DM at 95% C.L. These light parti¬ 
cles provide new decay channels for the SM Higgs boson, 
H}- Some of them contribute to the invisible branch¬ 
ing fraction, as is the case of the RH sneutrino, the 
RH neutrino and the lightest neutralino (H} —> NiNi, 
H} -a NN and H} —>■ XiXi)> whereas others would only 


HR {H} -a inv) 

BR (H° -A Non SM) 


= BR(iL 2 ° -A ViVi) 
+BR(iJ 2 ° -A NN) 
+BR (H} -> XiXi) 
= BR {H} -A H}H}) 
+BR {H} -a A}A}) 
+BR (H} -A inv). 


(4) 


Mostly only the RH sneutrinos contribute to the invisible 
branching fraction, since the RH neutrino mass is gener¬ 
ally too large. Neutralinos in this scenario would decay 
into the RH sneutrino through y) 1 — > Ni N or yj —> N\1 >l 
0 , contributing to the invisible branching ratio of the 
SM Higgs boson. Nevertheless, this contribution is neg¬ 
ligible in most of the points of the parameter space since 
either the lightest neutralino is too heavy or the branch¬ 
ing ratio of the H} -A x}xi process is very small. 


In Figure [Cl we plot BR(iL 2 -A Non SM) as a function 
of BR(7J 2 -a inv) for all the points that fit the GCE at 
95% C.L. We denote with dashed and dot-dashed lines 
the expected sensitivity in the future high luminosity 
LHC run to these two observables, respectively. Many 
solutions found here would lead to observable signals, 
mainly in searches of non SM decays of H}. Interest¬ 
ingly, this is the case of most of the solutions to the GCE 
that survive our estimation of the Pass 8 constraints on 
dSphs, which have BR(iJ 2 -A Non SM) > 0.1. Further¬ 
more, the A}A 5 final states that fulfil the Pass 8 bounds 
from dSphs might be also probed at the LHC, through 
the direct production of A} and its subsequent decay into 
r leptons producing a multilepton signal EzHU- 


spectrum, the better sensitivity of the Fermi-LAT to low-energies 

implies that the optimal value of E 7 is typically achieved close 8 Recall that the invisible decay modes are also included in the 
to the lower-end of the box. non-standard branching ratio. 
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FIG. 7. Left: m A o as a function of m H o for the points fulfilling all the experimental constraints and fitting the GCE at 95% C.L. with 
the colour convention as in Fig. [2] Right: Mass ranges for the Higgs sector and the different supersymmetric particles corresponding to 
the viable points in the parameter space. 


DISCUSSION OF THE VIABLE PARAMETER 
SPACE 

As seen previously, some of the best fits to the GCE 
in this model require the presence of either a light CP- 
even or a CP-odd Higgs boson (see Table U). In fact, after 
the inclusion of the dSph Pass 8 bounds, only the regions 
with a very light singlet-like CP-odd Higgs boson survive. 

In Figure [7] (left panel), we represent the solutions that 
fulfil all the experimental constraints and fit the GCE at 
95% C.L. in the plane (m^o, m H o). In general, points 
with light CP-even Higgs bosons below the mass thresh¬ 
old m H o /2 are more difficult to obtain than those with a 
light CP-odd Higgt@. The mass of the pseudoscalar Higgs 
can even be smaller than 2 mb (thus favouring its de¬ 
cay into leptons), whereas the mass of the lightest scalar 
Higgs boson remains above 20 GeV. It is noteworthy that 
the solutions allowed by the estimated Pass 8 bounds 
from dSphs always contain a pseudoscalar Higgs below 
50 GeV, and some of them also include a scalar Higgs 
below 60 GeV. 

In Figure [7] (right panel), we show the mass ranges 
of the different supersymmetric particles, as well as the 
Higgs sector for the aforementioned solutions. In gen¬ 
eral, these solutions have a mass spectrum that contains 


several non SM particles below the frontier of 100 GeV. 
Namely, we can see that neutralinos in these scenarios 
can be as light as 20 GeV while in some cases also the RH 
neutrinos can reach a mass around 60 GeV. Notice that 
the masses of sleptons and squarks do not vary much, 
since the soft masses and trilinear parameters for these 
particles were fixed in our scan. 

Figure [5] displays the solutions fitting the GCE at 95% 
C.L. in the planes k versus A (left panel) and tan/3 ver¬ 
sus /i (right panel), where the main constraint on these 
parameters is the bound on the Ul) mass [H]. For 
final states (blue) we observe a slight preference for small 
values of A, while for (cyan) and bb (grey) final 

states the solutions can be found for a wide range of val¬ 
ues. Very light pseudoscalar Higgses can be found in the 
Peccei-Quinn (PQ) limit n —> 0 and A K —> 0, for which 
the U(1)pq symmetry is restored, and it is precisely in 
these regions where we find most of the points that are in 
agreement with the expected dSph bounds using Pass 8 
data. In the right panel of Fig. [51 we can observe that 
these solutions are grouped in the region of \i < 180 GeV 
while the values of tan /3 do not show any preferred re¬ 
gion. The PQ limit is not needed for points that only 
fulfil Pass 7 data as the pseudoscalar mass can be larger. 


9 This is a consequence of the coupling of these Higgses to H § 
(the SM-like Higgs). Our choice of signs for A and re makes 
the coupling of H® to higher than the A® one; thereby the 
existing LHC constraints on the SM-like Higgs couplings to SM 
particles are more stringent for the H® final states I7fill . This does 
not apply to the limit re —> 0, where the relative sign between A 
and re does not play any role. 


Comparison with low-mass neutralino DM in the 
NMSSM 


The NMSSM can also accommodate low-mass neu¬ 
tralino DM 


77-8, 


Several studies have recently 
pointed out that these neutralinos can also explain the 
GCE 53-55J when they are either singlino/Higgsino or 
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FIG. 8. Left: k, versus A for the points that fulfil all the experimental constraints and fit the GCE at 95% C.L. with the colour convention 
as in Fig.[2] Right: The same but for the tan (3 and fi parameters. 



FIG. 9. m H o versus for the points fulfilling all the experi¬ 
mental constraints and fitting the GCE at 95% C.L. with the colour 
convention as in Fig. [2 The solid line corresponds to m H o = 
while the dashed line denotes the resonant annihilation condition, 
m H o = m N 1 / 2 ' 


bino/Higgsino admixtures [53| . The former situation can 
occur when k/ A « 1 and requires some degree of fine- 
tuning in the Higgs mass matrix parameters in order 
to avoid large deviations from the SM-like Higgs com¬ 
position. Moreover, to account for the observed relic 
abundance the mass of the pseudoscalar Higgs must be 
very close to the resonant condition, m^o ~ 2 m^o. The 
bino/Higgsino scenario requires k/X ^$> 1 and is however 
less fine-tuned since no resonant condition is required to 


fulfil the relic density constraint. 

The region of the parameter space in which RH sneu- 
trinos can explain the GCE is similar to that with 
singlino/Higgsino neutralinos, as k/X < 1 (in fact, light 
neutralinos are also frequent as we showed in Fig. 0 . It 
is worth noting that, especially when Pass 8 constraints 
are imposed, both models feature a similar Higgs spec¬ 
trum with light singlet-like scalar and pseudoscalar Higgs 
bosons. Nevertheless, in the NMSSM with RH sneutri- 
nos the amount of fine-tuning required is largely reduced 
as resonant annihilation is no longer necessary in order 
to obtain the correct relic abundance. In this model, the 
annihilation cross section increases when the annihilation 
channels into light scalar and pseudoscalar Higgs bosons 
(ViiVi —> H°H^ and iViiVi -> A°A °) are open [56]. This 
is illustrated in Figure 0 where m H o is represented as a 
function of m jy for the set of points for which the GCE 
is fit at 95% C.L. The solid line denotes m H o = m J y i 
(threshold for the annihilation in H/H/) and the dashed 
line corresponds to m H o = m^/2 (condition for reso¬ 
nant annihilatiorJR) . Only the points with bb final states 
require resonant annihilation. 

It should finally be emphasized that, contrary to the 
neutralino case, the parameters that determine the RH 
sneutrino properties (the soft mass, , the trilinear 
term, Ax,,, and the coupling, Xn) do not affect the 
masses of the Higgs sector of the NMSSM. Therefore, 
it is much easier to obtain the correct relic abundance 
without violating the LHC bounds. 


10 We remind the reader that the RH sneutrino does not couple to 
the light pseudoscalar. 
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CONCLUSIONS 

In this article, we have demonstrated that RH sneu- 
trino dark matter in the NMSSM can account for the 
observed low-energy excess in the Fermi-LAT gamma- 
ray spectrum from the Galactic Centre. Since we are 
working with a complete theoretical framework, we have 
also explored the complementarity with other DM search 
strategies, such as indirect searches for gamma-ray lines, 
direct DM detection, and the implications for a future 
LHC run. 

More specifically, we have performed a scan in the pa¬ 
rameter space of the model, incorporating all the con¬ 
straints from the LHC as well as direct and indirect dark 
matter searches. For the latter, we have included the 
Fermi-LAT bounds on dSphs and gamma-ray spectral 
features. In addition, we have estimated the effect of the 
latest Fermi-LAT reprocessed data (Pass 8). We have 
computed the gamma-ray spectrum for each point in the 
parameter space of the model, and then we have per¬ 
formed a x 2 fit to the excess. It should be emphasised 
that, contrary to usual model-independent approaches, 
we have taken into account the full annihilation prod¬ 
ucts, without assuming pure annihilation channels. 

We have obtained good fits to the GCE for a wide 
range of the RH sneutrino mass, ~ 20 — 135 GeV, 
with annihilation cross section in the DM halo £ 2 (cn>)o ~ 
5 x 10- 27 - 8 x 10 _26 cm 3 /s at 95% C.L. There is a 
large variety of final states for the RH sneutrino anni¬ 
hilation. In general, we observe that the fit to the GCE 
is good when the RH sneutrino annihilates mainly into 
pairs of light singlet-like scalar or pseudoscalar Higgs 
bosons, or A) 1 A 3 , that subsequently decay in flight. 

These produce spectral features, such as gamma-ray lines 
and/or box-shaped emissions that improve the goodness 
of the fit at large energies. The best fit of our analysis 
(y 2 = 20.8) corresponds to a RH sneutrino with a mass 
of 64 GeV which annihilates preferentially into a pair of 
light singlet-like pseudoscalar Higgs bosons (with masses 
of order 60 GeV). Good fits are also obtained for a va¬ 
riety of mixed final states and have been summarised in 
TablelH 

Our estimation of the constraints from the new Pass 8 
data from Fermi-LAT rules out all the points that fit the 
GCE at 68% C.L. and leaves only a few of them within 
the 95% C.L. region. Interestingly, these correspond 
mainly to points with predominant annihilation into light 
pseudoscalar Higgs bosons. It is remarkable that a re¬ 
gion with rrif,^ « 25 GeV remains valid.In this region, 
the light pseudoscalar has a small mass, m^o < 2 mb, 
and therefore leads to a leptonic final state with 4r, 
NiN\ —y 2A 3 —> 4r, which differs from the usual result 
for 2 t. 

We have then studied the implications for direct detec¬ 
tion of the points that fit the GCE at 95% C.L. Because 


of the presence of light non SM particles in the mass 
spectrum of the model, and the fact that many of our so¬ 
lutions involve resonant annihilation through the lightest 
scalar Higgs, the prospects for the RH sneutrino scatter¬ 
ing cross section off quarks span several orders of magni¬ 
tude. It is remarkable that many of the points found here 
are within the reach of next-generation experiments such 
as SuperCDMS and LZ. Interestingly, this is the case of 
the few points which also survive our estimation of the 
Pass 8 constraints. 

We have also investigated the prospects for the detec¬ 
tion of these scenarios in future searches for gamma-ray 
spectral features, using the estimated future sensitivity 
of Fermi-LAT for this kind of studies. We have found 
that only a small fraction of the points that fit the GCE 
at 95% C.L. will be observable in this way, mainly those 
associated with resonant annihilation through the light¬ 
est or second-lightest scalar Higgs boson. Unfortunately, 
none of the scenarios allowed by our estimation of the 
Pass 8 constraints would be detectable by spectral fea¬ 
ture searches. 

Finally, we have computed the contribution of these 
scenarios to the invisible and non-SM branching ratios 
of the SM Higgs and compared them with the predicted 
reach of the future high-luminosity LHC run at 14 TeV. 
We have determined that a substantial amount of viable 
scenarios (many of which also satisfy the Pass 8 con¬ 
straint) can be probed in such a way. 

It is remarkable that many of the scenarios that can 
explain the GCE with RH sneutrino DM can be subject 
to a complementary detection through other methods in 
the near future, such as direct DM searches or collider 
tests of the Higgs sector. 
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